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Abstract 



A measurement is presented of W-boson production in PbPb collisions carried out 
at a nucleon-nucleon (NN) centre-of-mass energy y/s^ of 2.76 TeV at the LHC using 
the CMS detector. In data corresponding to an integrated luminosity of 7.3 ub^ 1 , 
the number of W — > jiv^ decays is extracted in the region of muon pseudorapidity 
\t]P\ < 2.1 and transverse momentum pj > 25GeV/c. Yields of muons found per 
unit of pseudorapidity correspond to (159 ± 10(stat.) ± 12( syst.)) x 10~ 8 W + and 
(154 ± 10(stat.) ± 12(syst.)) x 10 -8 W~ bosons per minimum-bias PbPb collision. The 
dependence of W production on the centrality of PbPb collisions is consistent with a 
scaling of the yield by the number of incoherent NN collisions. The yield of W bosons 
is also studied in a sample of pp interactions at \/s = 2.76 TeV corresponding to an 
integrated luminosity of 231 nb -1 . The individual W + and W~ yields in PbPb and 
pp collisions are found to agree, once the neutron and proton content in Pb nuclei 
is taken into account. Likewise, the difference observed in the dependence of the 
positive and negative muon production on pseudorapidity is consistent with next-to- 
leading-order perturbative QCD calculations. 
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1 Introduction 

The hot and dense matter produced in heavy-ion (AA) collisions can be studied in a variety of 
ways. One approach is to compare AA to proton-proton (pp) collisions as well as to collisions 
of protons or deuterons with nuclei. Another way is to compare yields of particles whose prop- 
erties are modified by the produced medium to those of unmodified reference particles in the 
same AA collisions. Direct photons play the reference role at the Relativistic Heavy Ion Col- 
lider (RHIC) [lj and, more recently also at the Large Hadron Collider (LHC) |2H3]|. However, 
their measurement is complicated by copious background from tt° and t] meson decays, and by 
the presence of photons produced in fragmentation processes of final-state partons that can be 
affected by the medium (U. At LHC energies, new and cleaner references such as weak bosons 
in their leptonic decay modes become available (5HZI- The ATLAS and CMS collaborations 
recently reported first observations of Z bosons in heavy-ion interactions, showing that their 
yields per nucleon-nucleon (NN) collision are essentially unmodified by the medium DUEL 

Weak-boson production is recognised as an important benchmark process at hadron colliders. 
Measurements at 7TeV centre-of-mass (CM) energy in pp collisions at the LHC IfToHTZ l and 
previously, at other hadron colliders (Tevatron H3H3, RHIC EDI EH and SppS (2211231) with 
various collision energies, are well described by calculations based on higher-order pertur- 
bative quantum chromodynamics (QCD) using recent parton distribution functions (PDF). In 
PbPb collisions, W-boson production can be affected by initial-state conditions |5ll24lT26 1 . such 
as the mix of protons and neutrons. Since the leading-order W-production processes ud — > W + 
and du — > W~ reflect mainly interactions that take place between valence quarks and sea anti- 
quarks, the individual W + and W~ rates are expected to be modified relative to pp collisions, 
but not their sum. This is often referred to as the isospin effect, as it stems from a different 
content of u and d quarks in the proton relative to lead nuclei. The PDF can also be modified in 
nuclei, as parton depletion (or shadowing) could change the yield of W bosons at the LHC by 
as much as 15% in certain regions of kinematics [26j. Precise measurements of W production 
in heavy-ion collisions can therefore constrain the nuclear PDF and, moreover, provide insight 
into the PDF for neutrons. 

The W — > lv\ decays are of particular interest, since the charged leptons (/) lose negligible 
energy in the produced medium, regardless of its nature (partonic or hadronic) or specific 
properties J6H7|. Since they are dominantly created from a left-handed valence quark and a 
right-handed sea antiquark, W bosons are mostly left-handed and emitted in the valence quark 
direction, thus towards non-zero rapidity. The W + decays to a left-handed neutrino and a 
right-handed positive lepton, which is thus boosted back towards midrapidity while the W~ 
decays to a left-handed negative lepton which is boosted towards higher rapidity. This fact 
creates a difference in /+ and l~ yields as a function of lepton pseudorapidity rj, defined as 
rj = — ln[tan(0/2)], 9 being the polar angle of a particle trajectory with respect to the direc- 
tion of the anticlockwise-circulating heavy-ion beam. This angular difference and the relative 
abundances of W + and W~ bosons produced in PbPb compared to pp collisions (isospin effect) 
manifests itself in a lepton charge asymmetry, defined as a difference in Z + and l~ contributions 
divided by their sum. The measurement of this asymmetry as a function of muon pseudorapid- 
ity is quite robust, as it is insensitive to many systematic uncertainties. The W bosons appear 
therefore to be well suited to probe the characteristics of the initial state of PbPb collisions at 
LHC energies. 

This Letter reports the observation of W-boson production in a minimum bias (MB) sample 
of Nmb = 55.7 x 10 6 events from PbPb collisions collected with the Compact Muon Solenoid 
(CMS) detector at the CM energy for colliding nucleon pairs of y/s^ = 2.76 TeV. This sample 
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corresponds to an integrated luminosity of (7.3 ± 0.3) }ib~ l . These data were recorded during 
the first PbPb LHC data taking period at the end of 2010. In addition, we present results of a 
comparison analysis of W production in pp interactions in data obtained at the same y/s^ for 
an integrated luminosity of (231 ± 14) nb -1 , which is of a similar size to the nucleon-nucleon 
equivalent luminosity of the PbPb data. 

The Letter is organised as follows: the CMS detector is briefly described in Section 2, followed 
by the description of the experimental methods used for online and offline data selection in 
the PbPb and pp collected samples of events. The Monte Carlo (MC) simulations and the 
acceptance and efficiency correction factors derived from them are described there as well. 
The results and their discussion, together with comparison with theoretical predictions, are 
presented in Section 3. Finally the conclusions of this study are summarised in Section 4. 

2 Experimental methods 

A detailed description of the CMS detector can be found in Ref. Ij27| . In brief, a silicon pixel 
and strip tracker is located within a superconducting solenoid of 6 m internal diameter that 
provides a magnetic field of 3.8 T. The tracker consists of 66 million pixel and 10 million strip- 
detector channels, used to measure charged-particle trajectories for \tj\ < 2.5. It provides a 
vertex resolution of ~15 /mi in the transverse plane. Located within the solenoid, but outside 
of the tracker, are a crystal electromagnetic calorimeter and a brass/ scintillator hadron calori- 
meter. Muons are measured within \t]\ < 2.4 in gaseous detector planes embedded in the steel 
return yoke of the magnet. A matching of outer muon trajectories to the tracks measured in 
the silicon tracker provides a transverse momentum (px) resolution between 1 and 2%, for pj 
values up to lOOGeV/c. In addition, CMS has extensive forward calorimetry, in particular, two 
steel/ quartz-fibre Cherenkov, forward hadron calorimeters (HF), on each side of the collision 
point, covering 2.9 < \tj\ < 5.2. 

The centrality of PbPb collisions reflects the geometric overlap (impact parameter) of the in- 
coming nuclei, and is related to the energy released in these collisions and the effective number 
of NN interactions. CMS defines the centrality of a PbPb collision through bins that correspond 
to fractions of the total hadronic inelastic cross section, as observed in the distribution of the 
sum of the energy deposited in the HF 1128 H29 1 . The five bins in centrality used in this analysis, 
ordered from the smallest to the largest energy deposited in the HF, range from the most pe- 
ripheral, 50-100%, 30-50%, 20-30%, 10-20%, to the most central, 0-10%, collisions. These bins 
can be related through a Glauber model [30J to the number of nucleon-nucleon collisions per 
event. 

In this analysis, W bosons are measured through their W — > ]iVu decays. Muons can be cleanly 
identified and reconstructed, despite the high-multiplicity environment of heavy-ion collisions, 
a fact that makes this channel particularly suitable for measuring W production. The muon 
charge and transverse momentum vector are evaluated from the curvature of the track in the 
silicon tracker. The neutrino is not detected, but a large imbalance in the vector sum of the 
transverse momenta of all charged particles measured in the tracker is used to signal its pres- 
ence. 

A sample of MB events is selected that have a reconstructed primary vertex based on at least 
two tracks, and an offline-determined coincidence of energy depositions in both HF calorime- 
ters, with at least three towers, each above 3GeV. These criteria reduce contributions from 
single-beam interactions with the environment (e.g. beam-gas and beam-halo collisions within 
the beam pipe), ultra-peripheral electromagnetic collisions and cosmic-ray muons. The accep- 
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tance of this selection corresponds to (97 ± 3)% of the hadronic PbPb inelastic cross section [28J. 

Events for this analysis are selected using the two-level trigger of CMS. At the first (hardware) 
level, one muon candidate with a pj of at least 3GeV/c is required in the muon detectors. 
At the software-based higher level, one reconstructed track with a more precisely determined 
px > 3GeV/c is again required in the muon detectors. For muons from W-boson decays, the 
single-muon trigger efficiency is estimated as (97.0 ± 2.3)%. 

Muon offline reconstruction has k>99% efficiency to find tracks when hits in the muon detectors 
are taken as seeds. These tracks (called stand-alone muons) are matched to tracks reconstructed 
in the silicon tracker by means of an algorithm optimised for the heavy-ion environment 11291 
|31~] . For a muon from W decays, the silicon-tracking efficiency is ~85%, which is less than for pp 
collisions, as track reconstruction in the PbPb environment requires more pixel hits to reduce 
the number of possible combinations resulting from large particle multiplicities. Combined fits 
of the stand-alone muon and tracker trajectories (called global muons) are used in extracting 
the results of this analysis. Muon pseudorapidities are restricted to \rj^\ < 2.1, which provides 
uniform and good resolution both at the trigger stage and in offline reconstruction. 

A Z-boson veto is applied to reject events that contain a second muon of opposite charge with 
pj > lOGeV/c that forms a dimuon invariant mass of 60 < ntup < 120GeV/c 2 . Background 
muons from cosmic rays and heavy-quark semileptonic decays are rejected by requiring a 
transverse impact parameter of less than 0.3 mm relative to the measured vertex. No muon iso- 
lation criteria are required. The single-muon pj spectrum following this selection is shown in 
Fig. [T]ja) with red-filled circles. The enhancement in the number of muons with pj > 25 GeV/ c, 
expected from the decay of W bosons (green-hatched histogram), is evident. Details on the fit 
to the data are given below. 

To further characterise events with muons arising from W decays, the imbalance (p'j ) in the 
sum of the charged-particle transverse momenta with pj > 3GeV/c is computed for each 
event. The mean value of this transverse-momentum imbalance as a function of centrality 
of the PbPb collision is presented in Fig. [TJb) for data (black-filled squares) selected with the 
two-level muon trigger described above. The presence of significant pj in central events is ex- 
pected as these events contain many particles that are not included in the sum, such as neutrals 
or charged particles produced at low transverse momentum or at large pseudorapidity. For pe- 
ripheral collisions, the net fi tends to be quite small. Once a high-pj muon is required in the 
data (red-filled circles), the {p"j ) shifts to higher values of ~40 GeV/c, and is far less dependent 
on the centrality of the collision. This agrees with expectations (green triangles) for tfj values 
of undetected neutrinos originating from W decay. To enhance the contribution from the W 
signal, events are therefore required to have p l j > 25 GeV/c and pj > 20 GeV/c. 

The distribution in transverse mass (mj) for the p. and tfj system, computed as 

mj = J 2pjpx (1 — cos (p), where (p is the azimuth between the p\ and the p'j vectors, is pre- 
sented in Fig.JTJc) (red-filled circles), together with the expectation from a sample of simulated 
W events (green-hatched histogram) generated with PYTHIA v.6.4 [32] that are passed through a 
detailed simulation of the CMS detector based on Geant4 [33J. The effect of larger background 
from the underlying event in PbPb collisions is taken into account by embedding detector-level 
signals from the simulated W + and W~ decays into PbPb events generated with hydjet |34| 
(referred to as pythia+hydjet samples in the rest of the Letter). The mj spectrum of Fig.JTJc) 
is expected to have a sharp falling edge at the mass of the W boson, which is smeared by ex- 
perimental resolution. The region of mj > 40 GeV/ c 2 is used to define the W signal. Following 
this final selection, a total of 275 p + and 264 p~ events remain in the data sample. 
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Figure 1: (a) Single-muon transverse-momentum spectrum for |f/f | <2.1 in PbPb data (red- 
filled circles). Signal (green-hatched histogram) and background (blue-dashed histogram) con- 
tributions are fitted (black solid line) to the data, (b) Mean value of tfj for charged tracks as a 
function of centrality, before any event selection is applied on the muon-triggered data (black 
squares) and after it (red-filled circles), together with predictions from the PYTHIA+HYDJET 
samples (green triangles), (c) Transverse mass distribution for selected events in PbPb (red- 
filled circles) and pp (blue open squares) data, compared to simulation (green-hatched his- 
togram). The error bars represent statistical uncertainties. (See the text for more details). 
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Residual contributions from Z bosons with a misidentified muon or or a muon emitted in an 
insensitive region of the detector, as well as contributions from W — > tv t processes, where the 
t decays into }iv }l v T , are estimated using pp events simulated with pythia at the correspond- 
ing CM energy. A total estimated background contamination of 2.1% of the selected sample, 
based on the pp MC simulation, is subtracted from the data, as such electroweak processes 
are expected to scale with the number of elementary nucleon-nucleon collisions 1 30 1 . Remain- 
ing contaminations from a variety of QCD processes (mainly from semileptonic heavy-quark 
decays) in both pp and PbPb data, are estimated by extrapolating the nij distribution for both 
isolated and non-isolated, muon-enriched samples into the regions of signal. The estimate from 
pp data provides an upper limit on possible contamination of the PbPb sample, since parton 
energy loss (jet quenching) fl28| in heavy-ion collisions can only lower the yields relative to pp 
production. As a cross-check, the same method is applied to PbPb data. In both cases, upper 
limits on the contamination from QCD processes of 1% of the total selected PbPb sample are es- 
tablished, and no additional correction is applied for residual background, but a 1% systematic 
uncertainty is attributed to these sources of background. 

A cross-check on the contributions from signal and background to the selected sample is ob- 
tained by fitting the muon pj spectrum in Fig. [TJa) to two components: one arising from 
W — » fiVu signal (a template is taken from PYTHlA-simulated pp events), and another used 
for the background, modelled by a modified Rayleigh function [10J with 3 parameters. The 
best fit, shown in Fig. [lja), gives W boson yields that are in agreement to within 3% with the 
number of events found in the main analysis. 

Muon detector acceptance is evaluated with a sample of W events generated with pythia and 
CTEQ6L [35J PDF, that contains a weighted mixture of proton-proton, neutron-neutron and 
neutron-proton interactions, representing the nucleon content of Pb nuclei. 

Efficiencies for triggering, reconstructing, and selecting events are estimated using the PYTHIA+ 
hydjet samples previously discussed. Such embedded events, after being processed through 
full CMS trigger emulation and event reconstruction, reveal that track characteristics, such as 
the number of hits and the x 2 f° r the fit of the hits to muon trajectories, have similar distribu- 
tions in data and in simulation. The efficiency is evaluated separately for events with positive 
and negative muons. 

The acceptance parameter a is defined as the fraction of W bosons generated in the total phase 
space that decay into a muon with \r\V \ < 2.1 and pj > 25 GeV/c. For W + decays, a is estimated 
as 63%, compared to 49% for W~ events. The difference in the correction factors has its origin 
in the different angular distributions of the positive and negative muons from the W decays, 
as discussed in the Introduction. Within this acceptance, the overall trigger, reconstruction, 
and selection efficiency e, including the fj and nij criteria, averages to (73 ± 8)%, where the 
variations with the centrality of the event are considered in the uncertainty. To account for pos- 
sible detector effects, efficiency corrections are applied as a function of muon pseudorapidity 
separately for events with positive and negative muons. 

The individual components of muon efficiency are also estimated using data by means of a 
tag-and-probe method, discussed in Refs. HI HOfl , which entails counting Z candidates both 
with and without applying a probe requirement on one of the muons in Z -4 decays, 
to estimate: (i) the stand-alone muon-reconstruction efficiency, which is probed using charged 
tracks from the silicon tracker, (ii) the reconstruction efficiency of the silicon tracker, probed 
using stand-alone muons, and (iii) the trigger efficiency, probed by measuring the trigger re- 
sponse to the second muon in events triggered by a single-muon requirement. The latter is also 
checked with high-quality reconstructed muons found in minimum-bias events. In all cases, 
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these efficiencies, estimated with data, agree within the statistical uncertainties with those ob- 
tained from simulation. 

The total systematic uncertainty on the yield of W bosons is estimated as 7.4%, and is defined 
by summing separate contributions in quadrature as follows. The largest uncertainty is associ- 
ated with tracking efficiency, and corresponds to 4.9% accuracy obtained for the tag-and-probe 
determination of efficiency from data. Similarly, the uncertainty associated with the muon 
trigger is estimated as 2.3%. As indicated previously, the 1% maximum contribution from un- 
subtracted background sources is taken as a systematic uncertainty. The procedure used to 
estimate p"j has a 4% uncertainty, obtained by examining the impact of changing the threshold, 
from 1 to lOGeV/c, on the pj of tracks used to calculate the pj imbalance. The contribution 
from uncertainties on the calibration and resolution of the pj of muons is 0.2%. The probabil- 
ity of misidentifying muon charge is negligible (~10~ 4 ) and is ignored. The trigger efficiency 
for minimum-bias events is known to 3% accuracy. The measurement of W production is per- 
formed in the region of phase space defined by < 2.1 and p\ > 25GeV/c, but when the 
results are extrapolated to a larger region of phase space, using the acceptance corrections de- 
tailed above, the total systematic uncertainty rises to 8.2%. This is due to the 3.5% precision on 
the correction for acceptance that arises from uncertainties on the parameters of the assumed 
model (pythia) and the choice of PDF (CTEQ6L) used to create the weighted mixture of NN 
interactions, as well as uncertainties from higher-order QCD effects, ignored higher-order elec- 
troweak corrections, and the changes introduced in the underlying kinematics (y pj) of muons 
with the modifications made in the associated parameters. Table [l] presents a summary of the 
different sources contributing to the total systematic uncertainty in the analysis of PbPb data. 

Table 1: Sources of systematic uncertainties in the analyses of PbPb and pp data. The last line 
of the table contains the uncertainty in luminosity for pp data, and the analogous uncertainty 
in the value of Taa for MB events in PbPb data, as described in Section 3. 



Sources 


PbPb (%) 


PP (%) 


MB trigger 


3.0 




Muon trigger 


2.3 


2.0 


Tracking efficiency 


4.9 




Pj calibration 


0.2 


0.2 


Isolation 




1.3 


h 


4.0 


2.0 


Background 


1.0 


1.0 


Total experimental 


7.4 


3.5 


Acceptance 


3.5 


2.8 


Lumi (or equivalent) 


5.7 


6.0 



The analysis of pp events follows a similar procedure to that described for the PbPb data. Trig- 
gered events require the presence of a muon with a minimal value of pj, and those containing a 
muon reconstructed offline with pj > 25 GeV/c and \r\V \ < 2.1 are accepted, based on the same 
identification and quality criteria as for PbPb data. In addition, muon isolation is implemented 
by requiring the scalar sum of the transverse energies deposited in the calorimeters and of the 
transverse momenta of tracks (excluding the muon candidate) around the axis of the muon, 
within a cone of radius AR = a/ (A<£) 2 + (A?/) 2 = 0.3, to be below 0.1 p^. 

Unlike the transverse momentum imbalance based on charged tracks, used in PbPb events, the 
pr in pp events is calculated using a particle-flow (PF) technique [36j that combines tracking 
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and calorimetric information, with the requirement on fi of >20 GeV/ c implemented for the 
signal region. The resulting nij distribution for the muon and tfj system is shown in Fig. [TJc) 
(open blue squares), where reasonable agreement is observed with the signal reconstructed in 
PbPb data. The apparently-better resolution for pp relative to PbPb events can be attributed to 
the more powerful PF technique for measuring the p"j variable and, in general, to fewer parti- 
cles from the underlying event. After applying the same criteria as used in the PbPb analysis, 
a total of 301 }i + and 165 ]i~ events remain in the pp data. 

As in the PbPb analysis, the acceptance (a) and efficiency (e) for the inclusive pp — > WX — > 
jiv^X processes are evaluated using a sample of reconstructed W events in pp collisions at 
\Js = 2.76 TeV, simulated with pythia and processed through the CMS detector simulation. 
The mean values for acceptance are a = 61% for W + and tx = 54% for W~ events, and refer 
to W bosons produced in the entire phase space for W — > pv }l decays, with < 2.1 and 
Pj > 25 GeV/c. The average efficiency for selection within this acceptance is £ = 89%. 

The 3% background estimated in the selected pp sample is based on the same methods used in 
the PbPb analysis, and, as in the treatment of PbPb collisions, the contribution from electroweak 
background (2.1 %) is also subtracted frompp data. 

The systematic uncertainty on the evaluation of the pp cross section is estimated to be 7.0%, 
and, as before, obtained by summing in quadrature the uncertainty on luminosity (6%) and on 
the efficiency-corrected yield (3.5%), which is affected by sources similar to those discussed in 
the PbPb analysis: muon identification and trigger efficiencies (2%), isolation efficiency (1.3%), 
calibration and resolution of muon pj (0.2%), residual background (1%), and the calculation of 
pj (2%). Uncertainties on a correspond to 2.8%. Table [l] presents a summary of the different 
sources contributing to the total systematic uncertainty. 

3 Results 

Unless stated otherwise, all results reported in this section are evaluated in the restricted region 
of phase space defined by |//''| < 2.1 and p T > 25GeV/c. 

The yield (Nw) of muons from W decays per MB event in PbPb collisions, per unit of muon 
pseudorapidity is defined as Nw/Af/ = (N£ — N<P W )/ (cNmb^v)' where N^ el is the number 
of selected events, Ng W is the estimated background from other electroweak processes, Nmb = 
55.7 x 10 6 is the number of MB events, corrected for trigger efficiency, and estimated with an 
accuracy of 3%, e is the overall efficiency for W events, and At] = 4.2 is the pseudorapidity 
range used in the analysis. This yields N w +/A?/ = (159 ± 10 ± 12) x 10~ 8 and N w -/At] = 
(154 ± 10 ± 12) x 10~ 8 , for W + and W~ events, respectively, where the first uncertainty is 
statistical and the second is systematic. Extrapolating these measurements to the total W phase 
space using the acceptance corrections reported above, provides the total yields per MB event 
of N(W) = (1057 ± 63 ± 88) x HT 8 for W+, and (1317 ± 80 ± 108) x 10~ 8 for W~. 

The above analysis is repeated after subdividing the data into the five bins of event central- 
ity (defined in Section 2) and six bins in muon | rj \ . The total systematic uncertainty does not 
depend significantly on these variables, and is considered to be constant and uncorrelated be- 
tween bins. 

The yields of muons from W decays per MB event and per unit of muon pseudorapidity, 
Nw/ Atj, in PbPb collisions can be turned into the inclusive W cross sections per unit of muon 
pseudorapidity, normalised to the number of binary collisions occurring in PbPb interactions, 
when divided by the scaling factor Taa/ that is, (1/Taa) (Nw/Af/). This factor represents the 



8 



3 Results 



nuclear-overlap function, namely the number of elementary binary NN collisions divided by 
the elementary NN cross section, and can be interpreted as the NN-equivalent integrated lu- 
minosity per AA collision for a given centrality. In units of collisions per mb, the average Taa 
corresponds to 0.47 ± 0.07, 3.9 ± 0.4, 8.8 ± 0.6, 14.5 ± 0.8, and 23.2 ± 1.0, for event centralities of 
50-100%, 30-50%, 20-30%, 10-20% and 0-10%, respectively, and to 5.66 ± 0.32 for MB events, 
as computed using a Glauber model [30], with the parameters given in Ref. [28 1. The quoted 
uncertainties are obtained by changing the parameters of the model and the MB trigger and 
selection efficiencies by their respective uncertainties. 

Figure [2] shows the centrality dependence of the separate W + (violet-filled squares) and W~ 
(green-filled stars) normalised (NN-equivalent) cross sections and their sum (red-filled circles). 
The abscissa represents the average number of participating nucleons that undergo inelastic 
hadronic interactions (N par t) for the selected centrality intervals and computed using the same 
Glauber model. The open symbols at N par t ~ 120 correspond to MB events. For clarity, both 
W + and W~ points are slightly shifted relative to each other on the horizontal axis. Within 
present uncertainties, the normalised cross sections of W bosons in terms of the number of 
elementary nucleon-nucleon collisions are consistent with being independent of centrality of 
the PbPb collision, as can be seen by the results shown in Fig |2] and Table |2j 

Table 2: The selected number of events with }i + (N®+) and \C (N^ el ) and the normalised cross 

sections (1/Taa)(Mv/A//) as a function of centrality in PbPb data. The last row provides the 

selected number of events with ]d + (N f s ( + pp ) and }f~ (N^ el pp ) in pp data, and the cross sections 

derived from them (see text), divided by Arj. All values are reported for \t]V\ < 2.1 and pj > 
25GeV/c. 



Centrality (%) 


N sel 


Normalised cross section (W + ) 


N sel 

F 


Normalised cross section (W ) 






[nb]/A?/ 


[nb]/At] 


(50-100) 


11 


0.24 ± 0.07 ± 0.04 


10 


0.22 ± 0.07 ± 0.04 


(30-50) 


45 


0.30 ± 0.04 ± 0.04 


43 


0.28 ± 0.04 ± 0.04 


(20-30) 


48 


0.32 ± 0.05 ± 0.03 


46 


0.30 ± 0.04 ± 0.03 


(10-20) 


66 


0.27 ± 0.03 ± 0.02 


72 


0.28 ± 0.03 ± 0.02 


(0-10) 


105 


0.27 ± 0.03 ± 0.02 


93 


0.26 ± 0.03 ± 0.02 


(0-100) 


275 


0.28 ± 0.02 ± 0.02 


264 


0.27 ± 0.02 ± 0.02 


PP 


301 


0.34 ± 0.02 ± 0.02 


165 


0.18 ± 0.01 ± 0.01 



The W-production cross sections for pp collisions at the same \J s are determined in a similar 
manner, according to c pp = cr(pp — > WX) • B(W — > }iv pi ) = (Njf ~ pp — Ng W )/(e pp £), where 

B is the W boson leptonic decay branching fraction, Njf 1 pp is the number of selected events 
from the pp data sample, £ pp the overall efficiency for W bosons, and C the total integrated 
pp luminosity of 231nb _1 . The corresponding values, divided by the muon pseudorapidity 
interval of Arj = 4.2, are displayed at N part = 2 in Fig. [2] as open diamonds, and with the same 
colour code as for the PbPb data. The numerical values are given in Table [2] The impact of 
the neutron content in Pb nuclei is observed in the enhancement of W~ and the reduction of 
W + production in PbPb relative to pp interactions. As can be seen, the individual W + and W~ 
cross sections differ, while their sum agrees for PbPb and pp collisions. 

Theoretical predictions for W-boson production in PbPb collisions are based on the next-to- 
leading-order (NLO) MCFM |37l calculation and MSTW2008 PDF EH at next-to-NLO (NNLO), 
interfaced with the EPS09 [26J nuclear PDF that account for nuclear modifications in collisions 
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Figure 2: Centrality dependence of normalised W— > yiv^ cross sections (1/Taa)(Nw/At/) in 
PbPb collisions, for all W candidates (red-filled points) and, separated by charge, W + (violet- 
filled squares) and W~ (green-filled stars). The open symbols at N par t ~ 120 represent the MB 
events. At N part = 2, the corresponding cross sections are displayed for pp collisions divided by 
A77, for the same \fs. For clarity, both W + and W~ points are slightly shifted on the horizontal 
axis. The cross sections are given for the phase space region \ < 2.1 and pj > 25 GeV/c. The 
error bars represent the statistical, and the horizontal lines the systematic, uncertainties. 
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3 Results 



involving heavy ions rather than a sum over free nucleons. The effect of this nuclear PDF pack- 
age is a 4% reduction on the total W cross section with respect to the free nucleon PDF. The 
cross sections per NN pair, obtained in the region of phase space studied in this analysis, are 
(0.97 ± 0.10) nb for W+ and (0.87 ± 0.09) nb for W" production, where the quoted uncertain- 
ties take account of the two choices of PDF [26 1. Dividing these values by Arj = 4.2, provides 
predictions of (0.23 ± 0.02) nb for W + and (0.21 ± 0.02) nb for W~ production per unit of muon 
pseudorapidity values that are compatible with the experimental results given in Table [2] 

The cross sections for pp— > WX — > ]iv^X given above are now presented in the region of muon 
acceptance, denoted by <r acc ■ B, and are as follows: 



Pacc(pp -> w+x ) • B ( w+ -> H +v m) = (!- 44 ± °- 08 ± 0.10) nb; 
£7acc(pp^W-X)-B(W- ^}i-v F ) = (0.77 ± 0.06 ± 0.05) nb; 
^acc(pp -> WX) • B(W -> ftVjt) = (2.22 ± 0.10 ± 0.16) nb. 

Predictions for W-boson production in pp collisions at y/s = 2.76 GeV are obtained using the 
NNLO program FEWZ |3l| and the NNLO MSTW2008 PDF, and correspond to cross sections 
of (1.32 ± 0.06) nb for W + and (0.72 ± 0.03) nb for W~ in the same region of study. The uncer- 
tainties are from the choice of PDF. 

The experimental results extrapolated to the entire phase space of W-boson production, using 
the previous acceptance corrections, and referred to as cr tot • B, are: 



^tot(pp^W+X)-B(W+^+v F ) = (2.38 ± 0.14 ± 0.18) nb; 
0fet(pp -> W-X) • B(W~ -> p-Vp) = (1.45 ± 0.11 ± 0.11) nb; 
0tot(pp -> wx ) • B(W -> fiVp) = (3.83 ± 0.18 ± 0.29) nb. 

The corresponding NNLO predictions from FEWZ are (2.11 ± 0.10) nb for W+ and (1.29 ± 
0.05) nb for W~ production. The experimental results are compatible with the predictions, 
thereby confirming the validity of the standard model for W production at y/s = 2.76 TeV. 

The nuclear modification factors Raa = Nyy/ (Taa ■ Cpp)/ relating the W production in PbPb 
and in pp collisions, are computed from the measured yields in PbPb (Nw), the pp — > WX — > 
jiv^X measured cross sections, both quantities in the region of muon acceptance, and the 
nuclear-overlap function (Taa)- Although the overall W-production cross section is found to 
scale with the number of elementary collisions, the individual W + and W~ yields show a strong 
modification due to the nucleon content in Pb nuclei, as indicated by the extracted Raa factors 
for the region of phase space studied: 



Raa(W+) = 0.82 ± 0.07 ± 0.09; 
Kaa(W-) = 1.46 ± 0.14 ± 0.16; 
Raa(W) = 1.04 ± 0.07 ± 0.12. 

No cancellation of systematic uncertainties is assumed in computing these ratios. 

The difference in W + and W~ production at LHC and their subsequent leptonic decays pro- 
vide a different yield of W + — > }i + v^ and W~ — > defined as the "charge asymmetry" 
and given by A = (N w + — Nw- ) / (Nw+ + Nw- ), where Nw represents the efficiency-corrected 
number of selected events with a muon from a W decay, once the background from other elec- 
tro weak processes has been subtracted. Figure |3]shows this difference as a function of the muon 
pseudorapidity for PbPb collisions at 2.76 TeV (red-filled circles), for the experimental region 
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studied. The dot-dashed horizontal line at zero asymmetry is drawn for reference. Although 
the uncertainties are quite large, the measured asymmetry changes from positive to negative 
values, which indicates an excess of W~ over W + production at large \tj^\. The dependence of 
the asymmetry on muon pseudorapidity is in agreement with the predictions from MCFM and 
the MSTW2008 PDF, with small additional nuclear effects provided through the EPS09 nuclear 
PDF [26 1, and represented by the dashed curve. The uncertainty on the prediction associated to 
the use of both PDF ranges from ?»3% to ?»8% between small and large muon pseudorapidities. 

Most of the systematic uncertainties on these measurements affect ji + and }i~ events equally, 
and tend to cancel in the ratio A. A residual effect remains from a statistics-limited difference 
in efficiency observed for positive and negative muons in certain regions of pseudorapidity. 
The maximum difference (0.5%) is assigned as a systematic uncertainty. The effect of a possible 
difference in the calibration of muon pj or in the resolution for oppositely charged leptons 
is evaluated to be below 0.2%, and the impact of the subtraction of backgrounds from QCD 
processes is estimated to be of 1%. These factors amount to a common systematic uncertainty 
of 1.1% that affects each measured point. 

Results for pp collisions at 2.76 TeV are also shown in Fig. [3] (blue open squares), together with 
predictions from MCFM using the MSTW2008 PDF, represented by the solid curve. The uncer- 
tainty on the theoretical prediction is ~5%, which reflects the uncertainty from the choice of 
PDF. Central values obtained with MCFM, but using other PDF (CT10, CTEQ6.6M) gO), differ 
by !%10% from the result shown in Fig. [3| 

Effects due to the nucleon content in Pb nuclei are clearly visible in comparing results from 
PbPb and pp collisions. For the latter, the yield of }i + exceeds that of yT at all pseudora- 
pidities, reflecting the dominance of W + over W~ production. The integrated charge asym- 
metry A = 0.30 ± 0.04 ± 0.01 is in agreement with the MCFM prediction using MSTW2008 
PDF, which yields A = 0.30 ± 0.03. Recent measurements of this quantity in pp collisions at 
Vs = 7 TeV (12|, for the fiducial region analysed in this Letter, yield A = 0.189 ± 0.002 ± 0.008, 
indicating a dependence of the asymmetry on the CM energy of the interaction, which is in 
agreement with expectations from W kinematics and PDF evolution [38|. 

4 Conclusions 

The inclusive production of W bosons has been measured for W — > p.Vp decays in PbPb and 
pp collisions at ^/s^ = 2.76 TeV. Decays of W bosons were identified by requiring a recon- 
structed muon with \t]1*\ < 2.1 and p\ > 25GeV/c and transverse mass nij > 40GeV/c 2 . The 
W yields for all PbPb collision centralities were found to be consistent with those measured in 
pp collisions scaled by the corresponding number of incoherent nucleon-nucleon interactions. 
The individual W + and W~ boson yields are modified in PbPb compared to pp collisions due 
to the different proton and neutron content in the nuclear beams. The differences in the ji + 
and yields from W decays (charge asymmetries) have been measured as a function of muon 
pseudorapidity both in PbPb and pp collisions. All measurements were found to be well re- 
produced by higher-order perturbative QCD predictions. The pp results combined with those 
obtained at \Js = 7 TeV, exhibit the expected dependence of the charge asymmetry on the 
parton densities as probed at different collision energies. 

The results confirm the theoretical expectation that, in the probed ranges of parton fractional 
momentum and energy scale, further modifications of the nuclear parton distribution functions 
in the lead nucleus compared to the proton are small relative to the dominant isospin effect. 
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Figure 3: Charge asymmetry (Nw+ — Nyf-)/ (Nw+ + Nyf-) as a function of muon pseudora- 
pidity for PbPb (red-filled circles) and pp (blue open squares) collisions at ^/s = 2.76 TeV. 
Overlaid are predictions for pp collisions from MCFM calculations with MSTW2008 PDF (blue 
solid curve), as well as expectations for PbPb collisions from MCFM with MSTW2008 and EPS09 
PDF, that include isospin and nuclear effects (red dashed curve). Each prediction has an addi- 
tional uncertainty of 5%, estimated by the uncertainty from the choice of PDF. The experimental 
points have an additional 1.1% systematic uncertainty that is not shown in the figure. The dot- 
dashed horizontal line is drawn only for reference. 
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These studies demonstrate the promise of W-boson measurements as powerful tools in the 
investigation of initial and final-state effects in nuclear collisions at the LHC. The charge asym- 
metry of W-boson yields in PbPb interactions provides unique sensitivity to the parton dis- 
tribution functions for neutrons. Future analyses of larger data samples will yield enhanced 
constraints on the parton densities in nuclei and allow studies of W production in association 
with jets. 
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